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First observation of the decay 



The decay B~£ — > J / 'ipir + tt~ tt + is observed for the first time, using 0.8 fb _1 of 
pp collisions at yfs = 7 TeV collected by the LHCb experiment. The ratio of 
branching fractions B{B^ — > J/ipn^n ^tt + )/B{B^ — > J/ipTT + ) is measured to be 
2.41±0.30±0.33, where the first uncertainty is statistical and the second systematic. 
The result is in agreement with theoretical predictions. 
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The meson is the ground state of the be quark pair system^] Studies of its properties 
are important, since it is the only meson consisting of two different heavy quarks. It is also 
the only meson in which decays of both constituents compete with each other. Numerous 
predictions for B+ branching fractions have been published (for a review see e.g. Ref. (I)). 
To date, no measurements exist which would allow to test these predictions, even in ratios. 
.B+ production rates are about three orders of magnitude smaller at high energy colliders 
than for the other B mesons composed of a h quark and a light quark (B + , B° and 
Bg). On the experimental side, whatever is known about the B~£ meson was measured 
at the Tevatron. It was discovered by the CDF experiment in the semileptonic decay, 
B+ — > J/ipl + vX |2|. This decay mode was later used to measure the B+ lifetime [3j[4], 
which is a factor of three shorter than for the other B mesons as both b and c quark may 
decay. Only one hadronic decay mode of B+ has been observed so far, £> + — > J/ipn + . It 
was utilized by CDF |5) and D0 |6) to measure the B+ massQ 6277 ± 6 MeV 0. 

In this Letter, the first observation of the decay mode B+ — > J/ipiT + 7i^7i + is presented 
using a data sample corresponding to an integrated luminosity of 0.8 fb _1 collected in 
2011 by the LHCb detector |8j, in pp collisions at the LHC at y's = 7 TeV. The branching 
fraction for this decay is expected to be 1.5 — 2.3 times higher than for B+ — > J/iI>tt + (910 



However, the larger number of pions in the final state results in a smaller total detection 
efficiency due to the limited detector acceptance. We measure the B+ — > <//?/> 7r + 7r~7r + 
branching fraction relative to that for the B+ — > J/if)ir + decay and test the above theo- 
retical predictions. 

The LHCb detector [8] is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < r] < 5, designed for the study of particles containing b or c quarks. 
The detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4Tm, and three stations 
of silicon-strip detectors and straw drift-tubes placed downstream. The combined track- 
ing system has a momentum resolution Ap/p that varies from 0.4% at 5 GeV to 0.6% at 
100 GeV, and an impact parameter (IP) resolution of 20 [im for tracks with high transverse 
momentum. Charged hadrons are identified using two ring-imaging Cherenkov detectors. 
Photon, electron and hadron candidates are identified by a calorimeter system consist- 
ing of scintillating-pad and pre-shower detectors, an electromagnetic calorimeter and a 
hadronic calorimeter. Muons are identified by a muon system composed of alternating 
layers of iron and multiwire proportional chambers. The muon system, electromagnetic 
and hadron calorimeters provide the capability of first-level hardware triggering. The 
single and dimuon hardware triggers provide good efficiency for B^t — > J/ip7i + [7i^n + ], 
J/ip — > fi + fi~ events. Here, 7r + [7r~7r + ] stands for either tt + or 7r + 7r~7r + depending on 
the -B+ decay mode. Events passing the hardware trigger are read out and sent to an 
event-filter farm for further processing. Here, a software-based two-stage trigger reduces 



the rate from 1 MHz to about 3 kHz. The most efficient software triggers 11 for this 



analysis require a charged track with transverse momentum (px) of more than 1.7 GeV 



1 Charge-conjugate states are implied in this Letter. 
2 We use mass and momentum units in which c = 1. 



1 



(pr > 1-0 GeV if identified as muon) and with an IP to any primary pp- interaction vertex 
(PV) larger than 100 fim. A dimuon trigger requiring Pt(a<) > 0.5 GeV, large dimuon 
mass, M(yU + /i~) > 2.7 GeV, and with no IP requirement complements the single track 
triggers. At the final stage, we either require a J/0 — > fi + fi~ candidate with p^ > 2.7 GeV 
(> 1.5 GeV in the first 42% of data) or a muon-track pair with significant IP. 

In the subsequent offline analysis of the data, J/ip — > n + n~ candidates are selected 
with the following criteria: pr(fj) > 0.9 GeV, p T (J/0) > 3.0 GeV (> 1.5 GeV in 
the first 42% of data), y 2 per degree of freedom of the two muons forming a com- 
mon vertex, x^ tx (/i + yU~)/ndf < 9, and a mass window 3.04 < M(/i + /i~) < 3.14 GeV. 
We then find 7r + 7r~7r + combinations consistent with originating from a common vertex 
with x^ tx (7T + 7r~7r + )/ndf < 9, with each pion separated from all PVs by at least three 
standard deviations (XipC 71 ") > 9), and having Pt(tv) > 0.25 GeV. A loose kaon veto 
is applied using the particle identification system. A five-track J/?/>7r + 7r~7r + vertex is 
formed (Xvt x ( < W' 7r+7r ~7r + )/ndf < 9). To look for candidates in the normalization mode, 
B+ — > J/07r + , the criteria Pt(tt) > 1-5 GeV and Xvtx^/V' 71 " 4 ")/ 11 ^ < 16 are used. All 
candidates are required to have p^ > 4.0 GeV and a decay time of at least 0.25 ps. 
When more than one PV is reconstructed, that which gives the smallest IP significance 
for the B^t candidate is chosen. The invariant mass of a fi + {i~7t + [tt~7t + ] combination is 
evaluated after the muon pair is constrained to the J/ip mass and all final state particles 
are constrained to form a common vertex. 

Further background suppression is provided by an event selection based on a likeli- 
hood ratio. In the case of uncorrelated input variables this provides the most efficient 
discrimination between signal and background. The overall likelihood is a product of the 
probability density functions (PDFs), V(xi), for the four sensitive variables (xj): smallest 
Xip(vt) among the pion candidates, x^ x ( J/07r + [7r~7r + ])/ndf, -B+ candidate IP significance, 
xfp(B c ), and cosine of the largest opening angle between the J/ip and pion candidates in 
the plane transverse to the beam. The latter peaks at positive values for the signal as the 
.B+ meson has a high transverse momentum. Background events that combine particles 
from two different B mesons peak at negative values, whilst background events that in- 
clude random combinations of tracks are uniformly distributed. The signal PDFs, V S i g (xi), 
are obtained from a Monte Carlo simulation of B£ — > J/iIj7t + [tt~7t + ] decays. The back- 
ground PDFs, Pbk g ( obtained from the data with a <//?/> 7r + [7r 7r + ] invariant mass in 
the range 5.35 — 5.80 GeV or 6.80 — 8.50 GeV (far-sidebands). A logarithm of the ratio of 
the signal and background PDFs is formed: DLL sig / bkg = — 2 J2t=i ^(^sig^iV^bkgO^i))- 
Requirements on the log-likelihood ratio, DLL sig /bk g < —5 for £>+ — > J/ip n + 7Y^n + and 
DLL sig / bkg < —1 for B+ — > J/tpn + , have been chosen to maximize A^ sig / >jN sig + N hkg , 
where N sig is the expected B+ — > J/^7r + [7r~7r + ] signal yield and the Ab kg is the back- 
ground yield in the B+ peak region (±2.5 a). The absolute normalization of iV sig and A^kg 
is obtained from a fit to the <//?/> 7r + [7r~7r + ] invariant mass distribution with DLL S j g /bkg < 0, 
while their dependence on the DLL sig /bkg requirement comes from the signal simulation 
and the far-sidebands, respectively. The J/if)7c + [7T~7T + ] mass distributions after applying 
all requirements are shown in Fig. [T] To determine the signal yields, a Gaussian signal 
shape with position and width as free parameters is fitted to these distributions on top 
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Figure 1: Invariant mass distribution of —> J/ipir + TT tc + (top) and 
candidates. The maximum likelihood fits of B^ signals are superimposed. 



J/i/jtt + (bottom) 



of a background assumed to be an exponential function with a second order polynomial 
as argument. We observe 135 ± 14 B+ — > J/ip7i + 7i^7i + and 414 ± 25 B+ — > J/ipir + signal 
events. Using different signal and background parameterizations in the fits, the ratio of 
the signal yields changes by up to 3%. 

The ratio of event yields is converted into a measurement of the ratio of branching 
fractions B(B+ — > J/ip7i + 7i^7i + )/B(B^ — > J/ipir + ), where we rely on the simulation for the 
determination of the ratio of event selection efficiencies. The production of B+ mesons 
is simulated using the BCVEGPY generator |12, 13 1 which gives a good description of 
the observed transverse momentum and pseudorapidity distributions in our data. The 
simulation of the two-body B+ — > J/ip 7r + decay takes into account the spins of the particles 
and contains no ambiguities. The phenomenological model by Berezhnoy, Likhoded and 



Luchinsky 10 14] (BLL) is used to simulate B+ — > J/?/>7r + 7r~7r + decays. This model, 
which is based on amplitude factorisation into hadronic and weak currents, implements 
B+ — > J/ipW + * axial- vector form-factors and a W + * — > 7r + 7r~7r + decay via the exchange 
of the virtual af (1260) and p°(770) resonances. Since it is not possible to identify which of 
the same-sign pions originates from the p° decay, the two p° paths interfere. To explore the 
model dependence of the efficiency we also use two phase-space models, implementing the 
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Figure 2: Invariant mass distribution of the 7r + 7r _ 7r + combinations for the sideband-subtracted 
— > J/ip7r + Tr^ 7r + data (points) and signal simulation (lines). The solid blue line corresponds 
to the BLL simulations, the PH model is shown as a green dashed line and the PHPOL model 
is shown as a red dotted line. All error bars are statistical. 



same decay chain with no interference and with either no polarization in the decay (PH) 
or helicity amplitudes of 0.46, 0.87 and 0.20 for +1, and -1 J/0 helicities (PHPOL), 
respectively. For the helicity structure in the PHPOL model, we use the expectation 
for the B + — > _D*°a^(1260) decay based on QCD factorisation [15] . The background- 
subtracted distribution^] of the M(7r + 7r~7r + ) mass for the — > J/07r + 7r~7r + data shown 
in Fig. [2] exhibits an (2^(1260) peak and favours the BLL model. The p°(770) peak 
in the M(tt + ii~) mass distribution shown in Fig. [3] is smaller than in the two phase- 
space models, but more pronounced than in the BLL model, with the tail favouring the 
BLL model. The J/ip helicity angle distribution shown in Fig. [1] disfavours the model 
with no polarization. Since the BLL model gives the best overall description of the 
data, we choose it to evaluate the central value of the ratio of — > J/ip7i + 7i^7i + to 
B+ — > J/tp7i + efficiencies, 0.135 ± 0.004, and use the phase-space models to quantify the 

3 For comparisons between the data and simulation we use the data within ±2.5 a of the observed peak 
position in the mass (signal region) . We subtract the background distributions as estimated from the 
±(5 — 30) a near-sidebands. 
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Figure 3: Invariant mass distribution of the 7r + 7r~ combinations (two entries per candidate) 
for the sideband-subtracted — > J/ipTr + ir~iT + data (points) and signal simulation (lines). The 
solid blue line corresponds to the BLL simulations, the PH model is shown as a green dashed 
line and the PHPOL model is shown as a red dotted line. All error bars are statistical. 



systematic uncertainty. The phase-space models produce relative efficiencies different by 
—9% (PHPOL) and +5% (PH). We assign a 9% systematic uncertainty to the model 
dependence of — > J/V>7r + 7r~7r + efficiency. 

The distribution of the M(J/Y>7r + 7r~) mass has an isolated peak of four events at 
the ip(2S) mass. From the sidebands we expect 0.50 ± 0.25 background events in 
this peak. This is consistent with 3.6 ± 0.6 expected — > ip(2S)ir + events, assuming 
B(B+ ->■ iP(2S)tt+)/B(B+ J/V>7r+) equals to B(B+ ->■ ^(25)tt+) /B(B+ J/i/)ir+) = 
0.52 ± 0.07 1 7] after subtracting 10% to account for the phase-space difference. Since 
this contribution is only (2.6 ± 1.5)% of the B+ — > J/ip7i + ii^ii + signal yield, we do not 
subtract it and assign a 2% systematic uncertainty to the ratio of the branching fractions 
due to the efficiency difference between the B+ —> J/V , ai(1260) and B+ — > ip(2S)n + , 
ip{2S) — > J/ipir + TT~' decays, as obtained from the simulation. 

Other systematic uncertainties are due to limited knowledge of the B+ lifetime [7] 
(4%), uncertainty in the simulation of charged tracking efficiency (5%), trigger (4%) and 
the kaon veto (5%). Summing all contributions in quadrature, the total systematic error 
on the branching fractions ratio amounts to 14%. As a result, we measure the branching 
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Figure 4: Distributions of the cosine of the angle between the fi + and B + boosted to the rest 
frame of the J/tp meson for the sideband-subtracted B+ — > J/ifjir + (top) and B + — > J/tpTr + ir~ir + 
(bottom) data (points) and signal simulation (lines). In the bottom plot, the solid blue line 
corresponds to the BLL simulations, the PH model is shown as a green dashed line and the 
PHPOL model is shown as a red dotted line. All error bars are statistical. 



fraction ratio 

B(B+ M Ti+Ti- n + ) 

t / ^ = 2-41 ± 0.30 ± 0.33, 
B(B+ J/tp7i+) 

where the first uncertainty is statistical and the second systematic. 

The obtained result can be compared to theoretical predictions; these assume fac- 
torisation into B + — > J/ipW + * and W + * — > 7r + [/T~7r + ]. The contributions of strong 
interactions to B + — > J/ipW + * are included in form-factors which can be calculated in 
various approaches such as a non-relativistic quark model or sum rules. The coupling 
of a single pion to a W + * is described by the pion decay constant. The coupling of 
three pions to a W + * is measured in r~ — > u t 7v~tt + 7t~ decays, which are dominated by 
the ai(1260) resonance. The prediction by Rakitin and Koshkarev, using the no-recoil 
approximation in B+ J/0W+*, is B(B+ J/ipn + 7i-7i + )/B(B+ J/^tt+) = 1.5 
[9]. Likhoded and Luchinsky used three different approaches to predict the form fac- 
tors and obtained B(B + J li) tx+ tx~ tx+) j B{B+ ->- J/^vr+) = 1.9,2.0 and 2.3, re- 
spectively (To). Our result prefers the latter predictions. It is also consistent with 
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B{B + -»■ D*°7i + 7i-7i + )/B(B + -»■ D*°7r+) = 2.00 ± 0.25 0, which is mediated by similar 
decay mechanisms, and with a similiar ratio of phase-space factors. Our result constitutes 
the first test of theoretical predictions for branching fractions of B~£ decays. 
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